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Abstract 

The general problem of representing the Greens function G{k, z) in terms of self energy, in field 
theories lacking Wick's theorem, is considered. A simple construction shows that a Dyson like 
representation with a self energy S(A;,z) is always possible, provided we start with a spectral 
representation for G(/c, z) for finite systems and take the thermodynamic limit. The self energy 
itself can then be iteratively expressed in terms of another higher order self energy, capturing the 
spirit of Mori's formulation. 

We further discuss alternative and more general forms of G{k, z) that are possible. In particular, 
a recent theory by the author of extremely correlated Fermi liquids at density "n" , for Gutzwiller 
projected non canonical Fermi operators obtains a new form of the Greens function: 

z-Ek- ^{k, z) 

with a pair of self energies ^{z) and ^{z). Its relationship with the Dyson form is explored. A 
simple version of the two self energy model was shown recently to successfully fit several data sets 
of photoemission line shapes in cuprates. We provide details of the unusual spectral line shapes 

that arise in this model, with the characteristic skewed shape depending upon a single parameter. 
The EDO and MDC line shapes are shown to be skewed in opposite directions, and provide a 
testable prediction of the theory. 



I. INTRODUCTION 



Our work explores the representation of the Greens function G{k, z) of a particle, in 
field theories without Wick's theorem. While Wick's theorem for Bosons and Fermions 
automatically gives rise to the Dyson form of self energy, its absence for non canonical, i.e. 
general operators (other than Bosons of Fermions), leads to a conundrum that is poorly 
understood. This work addresses a particular type on non canonical theory originating from 
Gutzwiller projection of electrons on a lattice. 

The representation of the propagator or the Greens function in terms of its Dyson self en- 
ergy is a fundamental paradigm of standard interacting relativistic and non relativistic field 
theories. The structure of this representation and the generation of approximate approxima- 
tions for the self energy in terms of Feynman diagrams, the vertex function or higher order 
Greens functions form the dominant part of existing literature of many particle physics. 

In the context of extending these studies to extremely large and singular repulsive in- 
teractions, termed extreme correlations [T], one needs to deal with non canonical electrons. 
A standard non-canonical problem involves the Hubbard operators [2] X'^'^ located at sites 
j of a lattice. These are "graded" projection operators with label a representing the three 
allowed local configurations 0, t, 4- Among these operators, X^" , XJ^ are Fermionic destruc- 
tion and creation type objects. Their Greens function is measured directly in angle resolved 
photoemission experiments (ARPES) on certain experimental systems embodying extreme 
correlations, including the high temperature superconductors 0. 

Quite recently the author has formulated in Ref. (jTj), the theory of an Extremely Cor- 
related Fermi Liquid (ECFL) state of the t-J model, where he has found another type of 
representation for the Greens function with a pair of self energies (see Eq. ^ below) by 
using the non perturbative Schwinger approach of source fields to depict the equations of 
motion. The use of more general forms of Greens functions is not completely new, there 
are examples in literature of multiple self energies in Ref. (jl]) and Ref. (jSj). The physics 
of extreme correlations treated here, is based on non perturbative considerations without 
obvious parallels in weak or intermediate coupling problems. It leads to the two self energy 
form Eq. (|4]), whose distinctive signatures are strikingly different from those of Fermi liquids. 

The technical details of the construction in Ref. (jTj) are intricate and require the processing 
of two parallel hierarchies for the two self energies. A separate paper is in preparation 
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detailing the involved technical details and the calculation in Ref. ([T]). Further background 
details of the notation, definitions and sum rules satisfied by the Greens functions for extreme 
correlations, and its analyticity can be found in the earlier publication Ref. A suggestive 



functional form of the Greens function Eq. (37) emerges from Ref. (II]), by making the 
assumption of momentum independence of the two self energies, valid in high dimensions. It 
satisfies the number sum rule and the total particle weight integrates to unity in each state. 
In Ref. (j3]), Gweon, Gu and the author have shown that several experimental data sets on 
ARPES by different groups using both the traditional synchrotron source light and laser 



source light can be reconciled very well with the line shape in Eq. (37) and Eq. (60). This 
is the first satisfactory functional form that has been found to fit both laser and synchrotron 
data, and works very well with few adjustable parameters. 

Therefore a major objective of this paper is to elucidate the detailed form of the spec- 
tral lines that emerge from the above simple version of the ECFL theory in Ref. (fl]) and 
successfully employed to understand experiments in Ref. Our hope is that this de- 

tailed analysis would familiarize readers with the nuances of the new spectral function, and 
thereby facilitate ARPES line shape analysis of further experiments on high temperature 
superconductors and other materials, in a manner analogous to that in Ref. (j3]). 

For the above purpose, we recall that in a lattice of finite (say small) number of sites, the 
state space is finite dimensional and hence the Greens function for arbitrarily complicated 
objects can be computed by numerical means, leading to rational functions of the complex 
frequency z as in Eq. ([s]) below. We begin by studying this representation and see how 
the Dyson representation arises; we find that the two self energy representation Eq. ^ is 
also quite natural from this view point. We further study the infinite size limits where the 
poles coalesce to give cuts in the complex z plane. We adopt a phenomenological model for 
an underlying auxiliary Fermi Liquid (aux-FL) self energy, enabling us to display analytic 
expressions for the Greens function. We provide a detailed perspective on the representation 



in Eq. (37), namely the location of the poles and the subtle differences from a standard Fermi 
liquid. 

Another result in this paper is to show that a Dyson like representation with a self energy 
S(A;, z) is always possible, provided we start with a spectral representation. The self energy 
itself can then be iteratively expressed in terms of another higher order self energy. This 
hierarchical result is cast in the same form as the Mori formalism. While the Mori formalism 
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is very abstract, and expressed in terms of projection operators, we can go beyond it in a 
certain sense. By working with standard spectral representation, we show that it is possible 
to express the higher self energy spectral functions in terms of the lower ones, leading to 
an explicit hierarchy. Our construction completely bypasses the Mori projection operators, 
and should be useful in throwing light on the latter. 

The plan of the paper is as follows. In Sec (II) and (III) we note the spectral representa- 
tion and study the Greens function as a rational function of complex frequency z for a finite 
system. In Sec (IV) we note the representation in the limit of infinite size and introduce 
the high dimensional expression with two self energies. The detailed structure of the char- 



acteristic line shape as in Eq. (60) is discussed, and its dependence on physical parameters 
displayed with the help of a numerical example. An explicit example of the auxiliary Fermi 
liquid Greens function is provided, and typical values of the parameters are argued for. In 
Sec (V) the line shapes in EDO and MDC are displayed in detail, in order to bring out the 
specific signatures of the theory, namely a skew in the spectrum arising from the caparison 



factor in Ref. ([T]) and Eq. (37). In Sec(VI) we discuss the amusing connection with higher 
order self energies of the type that Mori's formalism yields, but at a much more explicit 
level than what is available in literature. 



II. SPECTRAL REPRESENTATION OF THE GREENS FUNCTION 

Let us begin with the spectral representation [10] of the Matsubara Greens function at 
finite temperatures given by: 

G{k,z)= [ dx^^^, (1) 
J z-x 

where G is the Greens function at a fixed wave vector fc, and pc is its spectral density and 
the integration range is — oo < x < oo. To simplify notation, we call the Greens function 
as G{k,z), the same object was denoted by Q{k^z) in Ref. (jT|. The index k can be also 
replaced by a spatial index when dealing with a local Greens function. The spectral function 
PG{k,x) in most problems of interest in Condensed Matter Physics has a compact support, 
so that G{k,z) has "reasonable" behaviour in the complex z plane, with an asymptotic 1/z 
fall off, and apart from a branch cut on a portion of the real line, is analytic. The frequency 
z is either Fermionic or Bosonic depending on the statistics of the underlying particles. The 
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spectral function is given by the standard formula[2l [TO] 



Pg{Kx) = ^\{a\A{k)\(3)\^{Pa.+ Pp) + (2) 

where A{k) is the destruction operator, is the Boltzmann probabihty of the state a given 
by e^^^" /Z and Ea is the eigenvalue of the grand Hamiltonian of the system K = H ~ fiN. 
In the case of canonical particles A{k) is the usual Fermi or Bose destruction operator. 
In Ref. (jl]), non canonical Hubbard "X" operators are considered; we will not require any 
detailed information about them here except that the anticommutator |y4, A"''} is not unity, 
but rather an object with a known expectation value {l — n/2), in terms of the dimensionless 
particle density n. 

We consider two alternate representations of the Greens function in terms of the complex 
frequency z that are available in many-body physics: (a) for canonical Bosons or Fermions 
the Dyson representation in terms of a single self energy and (b) for non canonical 

particles a novel form proposed recently by the author with two self energy type objects 
and "^{z). 

G{k,z) = ^ — , Dyson (3) 



z- Ek- Tj{k,z) 



ECFL (4) 



z- Ek- ^{k,z) 

For canonical objects = 1 and for Hubbard operators in the ECFL we write ac = l — n/2. 
We start below from a finite system, where the Greens function is a meromorphic function 
expressible as the sum over isolated poles in the complex frequency plane with given residues. 
In fact it is a rational function as well, expressible as the ratio of two polynomials. Using 
simple arguments, we will see that the above two representations in Eq. ^ and Eq. ^ are 
both natural ways of proceeding with the self energy concept. In the limit of a large system, 
the poles coalesce to give us cuts in the complex frequency plane with specific spectral 
densities. In this limit we display the equations relating the different spectral functions. 

III. FINITE SYSTEM GREENS FUNCTION 

We drop the explicit mention of the wave vector k, and start with the case of a finite-size 
system, where we may diagonalize the system exactly and assemble the Greens function from 
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the matrix elements of the operators A and the eigenenergies as in Eq. (|2|. We see that pc 
is a sum over say m delta functions located at the eigenenergies Ej (assumed distinct), so 
we can write the meromorphic representation 

m 

j=l 1 

The overbar in G{z) is to emphasize that we are dealing with the finite-size version of the 
Greens function G{z). Here CLj,Ej constitute 2m known real parameters. The sum 

m 

J2aj = aG, (6) 

i=i 

where ac = 1 for canonical objects and we denote = 1 — n/2 for the non canonical case 
of ECFL. In the infinite size limit we set G{z) — t- G{z). It is clear that for z ^ {Ej}max we 
get the asymptotic behaviour G — )■ and therefore G is a rational function that may be 
expressed as the ratio of two polynomials in z of degrees m — 1 and m: 



G{z)=aG^^, Q{z) = ll{z-E,), P{z)=ll{z-^r), (7) 

^ ' j=l r=l 

where the roots jr are expressible in terms of aj & Ej. We use the convention that all 
polynomials Qm ■ ■ ■ have the coefficient of the leading power of z as unity, and the degree is 
indicated explicitly. 

We now proceed to find the self energy type expansion for G, and for this purpose 
multiplying Eq. (|5]) by 2 and rearranging we get the "equation of motion" : 

{z-E)G{z) = aG + I{z), (8) 

where we introduced a mean energy E 



m—l 



m 



liz) = (9) 

j=l ^ 3 

SO that asymptotically at large z, we get I{z) ~ 0(l/z^). In standard theory E plays 
the role of the Hartree Fock self energy so that the remaining self energy vanishes at high 
frequencies . Motivated by the structure of the theory of extremely correlated Fermi 
systems Ref. m, we next introduce the basic decomposition: 

I{z) = G{z)<^{z) + ^!{z), (10) 



where we have introduced two self energy type functions $(-2) and "^{z) that will be deter- 
mined next. Clearly Eq. (10) leads immediately to the Greens function Eq. Q (or Eq. ([s]) 



if we set \1' — )■ 0). The rationale for Eq. (10) lies in the fact that the function / has the same 



poles as G{z). Thus it has a representation as a ratio of two polynomials 

with Rm-2 a polynomial of degree m — 2 and io a suitable constant and the same polynomial 
Q from Eq. ([7|, thereby it is natural to seek a proportionality with G itself. If we drop 
and rename $ — t- S, then this gives the usual Dyson self energy determined uniquely 
using Eq. ([T]) and Eq. (11) as: 



The expression Eq. (10) offers a more general possibility, where $(2) and "^{z) may be 
viewed as the quotient and remainder obtained by dividing I{z) by G{z). It is straight- 
forward to see that "^{z) and ^{z) are also rational functions expressible as ratios of two 
polynomials 

,T,/ \ I ^m-3{z) X , Lm-2{z) ,^ „s 

where K, L, D are polynomials of the displayed degree. Comparing the poles and the zeros 



of G beween Eq. \A\ with Eq. (12) and Eq. (^7|, we write down two equations 



Qm = {Z- E)Dm-l - (i)oLm-2, (14) 

SO that we may eliminate D and write an identity 

{Z - E)Pm^i - Qm = —{Z - E)Km-3 + (poLm-2- (15) 

Here the LHS is assumed known and we have two polynomials to determine from this 
equation. Therefore there are multiple solutions of this problem, and indeed setting K ^ 
gives the Dyson form special case. 



A. A simple example with two sites: 



The Greens function of the t-J model at density n with J = and only two sites is 
a trivial problem that illustrates the two possibihties discussed above. The two quantum 
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numbers k = 0,7i correspond to the bonding and antibonding states with energies = =F^, 
and a simple calculation at a given k gives Eq. ([s]) as 

Gik,z) = ^ + ^, (16) 
z -ek z + Ck 

where z = icj„ + /i, 02 = e^^{l + e^^^~'^''^) /{2Z), oi = 1 — n/2 — 02 and the grand partition 
function Z =1+ Ae^'^'' + Ae^^" cosh(/3t). This can be readily expressed as 

= z-Ek-m^) ' ^ ' = ^ ' = ^^^^ 

where Ek is arbitrary, = [E^ — el) and Bk = {1 — n/2){Ek — Ek) and with the first 
moment of energy Ek = ek{ai — a2) / {1 — ^1/2) . As we expected, the functions $ thus have 
a single pole, as opposed to G with two poles. In this case the dynamics it rather trivial so 
that the choice of Ek is free. If we set Ek = Ek the residue Bk vanishes and so the second 
form collapses. 



B. Summary of analysis: 

In summary, guided by analyticity and the pole structure of G{k,z), we find it possi- 
ble to go beyond the standard Dyson representation. However, we end up getting more 
freedom than we might have naively expected. This excess freedom is not unnatural, since 
we haven't yet discussed the microscopic origin of these two self energies. The theory in 
Ref. ^ provides an explicit expression for the two objects \1/ and $, where a common lin- 
ear functional differential operator L generates these self energies by acting upon different 
"seed" functions as in Eq. (7) of Ref. The above discussion therefore provides some 
intuitive understanding of the novel form of the Greens function in Eq. Q, without actually 
providing an alternative derivation to that in Ref. (jT]). 



IV. INFINITE SYSTEM SPECTRAL DENSITIES AND RELATIONSHIPS 

In the infinite size limit, the various functions will be represented in terms of spectral 
densities obtained from the coalescing of the poles. Following Eq. ([T]) we will denote a 
general function 

Qiz) = [ dx (18) 
z — X 



where Q = S, $, in terms of its density Pq{x). The density is given by Pq{x) = 
(— ^)53m Q{x + zO+), as usual. In parallel to the discussion of Eq. ([l]), the assumption 
of a compact support of pq gives us well behaved functions. We now turn to the objective 
of relating the spectral functions in the two representations discussed above. 



A. Spectral representation for the Dyson Self Energy 

Let us start with Eq. ([T]) and the standard Dyson form Eq. ([s]) where we drop the overbar 
and study the infinite system function G{z). We use the symbohc identity: 

^ V--i7iS{x), (19) 



X + zO+ X 

with real x, V denoting the principal value, and where the Hilbert transform of a function 
f{u) is defined by 

'H[f]i^)='P r dy^. (20) 

J -oc ^ y 

We note the following standard result for completeness: 

Pg{x) = ac (21) 

(vrpE(x))2+ (x-E-?/[ps](x)j 

A more interesting inverse problem is to solve for Pt,{x) given G{z). Towards this end we 
rewrite the Dyson equation as 

m = ^-E-^^, (22) 
where the self energy vanishes asymptotically as 1/^ provided the constant part, if any, is 



absorbed in E. Therefore this object can be decomposed in the fashion of Eq. (18). We 



compare Eq. (22) with Eq. (18) with Q — S and conclude that 



Pii[x) = — — = 7T. (23) 

vr G{x + tO+) {7ipG{x)y + {^eG{x)f ^ ' 

The real part can be found either by taking the Hilbert transform, 

^eT.{x) ='H[p^]{x), (24) 

or more directly as 

^enx)=x-E-^e—^^ = x-E ^_g3^^ (25) 

^(^ + ^0+) {'KpG{x)Y + {^eG{x)f 
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B. Spectral representation for the ECFL self energies: 

For the ECFL Greens function in Eq. (4), we set ac = (l — f) and write E ^ ^ 
representing the single particle energy measured from the chemical potential. We start 
with the expression: 

and express it in terms of the two spectral functions and p$ [13]. We can write spectral 
function pc'- 

where A(,^, x) and the term rj are defined as: 

i)(f,i) = Mbj.|(f,J^) + ^|i^WM«,i). (29) 

The real part of G is also easily found as 

. ^ ((1 - f) + n[p^u,x)) {x-i- n[p^u,x)) - 7rV^(e,x)p^(e,x) 

(vrp$(e,a;))2 + (x-^-H[p*](e,x)f 
Thus given the ECFL form of the Greens function, we can calculate the Dyson Schwinger 



form of self energy in a straightforward way using the inversion formula Eq. ( 23 ) and Eq. ( 24 ) . 
The inverse problem of finding $ and \E' from a given S or G is expected to be ill defined, 
as discussed above for finite systems. 



The first Fermi liquid factor in Eq. (|27|) has a peak at the Fermi liquid quasiparticle 

'F 
k 



frequency E^^ for a given given as the root of 



-ik- H[pMk, = 0, (31) 

however pa itself has a slight shift in the peak due to the linear x dependence in the 
numerator. This is analysed in detail in the next section for a model self energy. At this 



solution x(^), Eq. (30) gives a relation: 

p*(a, i^D = -P*(6, X 3fte (32) 
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C. High Dimensional ECFL Model with k independent self energies and its Dyson 
representation: 



In this section we illustrate the two self energies and their relationships in the context 
of the recent work on the ECFL Ref. ({T]), and in Ref. ([3]). Here we study a model Greens 
function, proposed in Ref. for the t-J model, that should be suitable in high enough 
dimensions. It is sufficiently simple so that most calculations can be done analytically. 
The model Greens function satisfies the Luttinger Ward sum rule[B] and thereby maintains 
the Fermi surface of the Fermi gas, but yields spectral functions that are qualitatively 
different from the Fermi liquid. This dichotomy is possible since it corresponds to a simple 
approximation within a formalism that is very far from the standard Dyson theory, as 
explained in the previous sections. Our aim in this section is to take this model Greens 
function of the ECFL and to express it in terms of the Dyson self energy so as to provide a 
greater feel for the model. 

Here the two self energies are taken to be frequency dependent but momentum indepen- 
dent, and using the formalism of Ref. become related through Aq, an important physical 
parameter of the theory: 

m = -iA^-^w- (33) 

The physical meaning of Aq as the mean inelasticity of the auxiliary Fermi liquid (aux-FL) 
is emphasized in Ref. ([T]), and follows from Eq. (56). Thus = —^^P<!>, and hence we 
get the simple result [TT] : 

,2 



1 \ 

2/ 4Ao ^ ^ 



(34) 



The auxiliary Fermi liquid has a Greens function g~^{C,k, z) = z — C,k ~ ^{z), where is the 
electronic energy at wave vector k measured from the chemical potential p, and therefore 
we may write the model Greens function as 

where 

4 , , . 

= Ao ^ (1 - -). (36) 

With r(x) = 7rp$(x), 3?e <l>(x + iO^) = 'H[p$](x) and e{iki = {x — C,k — 'H[p#](x)), we can 
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0.10 



FIG. 1: The density n = .85, temperature T = 600K, Aq = 0.0786 and parameters are Set(I) 



of Eq. (54). At this rather high temperature, we can see the details of the spectral shape clearly. 



The vertical line is at x* = E^^, this energy is the location of the peak of the physical spectral 
function pc as marked. Its leftward (i.e. red) shift relative to the Fermi liquid peak at the chemical 
potential is clearly seen. The two horizontal lines specify the magnitude of the 5Re G{0, x) at x = 
(HI) and x = E*^ (H2). The line HI is at height nV(4Ao) and H2 is at height nV(4Ao)(l -Zfe/2). 



express the spectral function and the real part of the Greens function as 

^2 ^ T(x) 



(47rAo; T\x)+e\ik,x) 



{£^0 + - x} 



n 

4Ar 



1 + 



£i^k,x){€o + -X) 

T^{x)+e^i^k,x) 



(37) 
(38) 



The linear frequency term in braces in Eq. (37) is termed the caparison factor in Ref. 



and leads to significant features of the spectrum as discussed below. For completeness we 
note the auxiliary Fermi liquid part of the problem as 



Psi^k,x) 



1 



(39) 



7TT^ix)+e''{^k,x) 

In Fig. ([T]) we plot the above three functions for a model system described more fully in 
Section( [lVF| . 



D. EDC or constant wave vector scans and energy dispersion. 
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We first study the peak structure corresponding to fixing k the wave vector and hence 
^k, and sweeping the energy x. These give rise to the energy distribution curves, i.e. the 
EDO's. The aux-FL part has a peak at x = E^^ for a given ^fc, as in standard FL theory 
from solving for the roots of Eq. (31 ). For A; ~ fci? we find E^^ = ^kZk with the momentum 



independent self energy ^{z), where 

Zk= lim {l-d^e(^{x)/dx)-\ (40) 
Expanding around this solution we write 

£{^k,x)^^{x-Zk^k), (41) 



We will also write F^. = T{x)/ ^_^^fl at the FL quasiparticle location, where we expect 

k 

for the Fermi liquid F^ ~ ci{k — kpY + C2T^, with suitable values as described more fully 



in Section( IV F). At this value, we have the identity 3fte G{^k,E^^) = as remarked 



above. As a consequence, in Fig. ([T]) the intersection of the line HI and the vertical y axis 
also coincides with the value of 3fte G at the chemical potential. To elucidate the line shape 
of the ECFL, we start with the FL solution and perturb around it to find the corrected 
location of the peaks in the full spectral function. 

Similarly the real part is found 

^2 r _ Ei^) {So + - xY 



iPeak I - - ^ 



zin + ix-Ei^Y 



4Ao 

We introduce the following convenient positive variable Uk- 



(43) 



■ _ ^0 + ik- 

smhwfc = — , (44) 

^fc-L k 

SO that near the Fermi energy and at low T the small F^ ~ drives it to a large and 
positive value i.e. exp-u^ — )• 1/T^ . At higher binding energies, decreases towards zero. 



as discussed below. We will also define a dimensionless variable Qk below in Eq. (61) that 
depends on Uk only and determines the shape of the peak. To analyze the shape at a given 
we introduce a dimensionless energy variable e through the relation 

X = ^ + ZkVk e, (45) 
13 



T=300K 



Peak ratio 




-0.05 



FIG. 2: The density n = .85, T = 300K, Aq = 0.0678 and parameters are from Set(I) in Eq. (54). 



The various dimensionless variables shown against E^^ are the peak ratio from Eq. (50), the 



variable Uk (scaled by 10) from Eq. (44), the skew asymmetry factor /t(i^fc) from Eq. (65) and the 
variable from Eq. (|61|). 



where we must require that |e| ~ 1 for the expansion around the FL peak to be valid. The 
spectral function is expressible as 

sinh(-Ufc) — e 



and 



n 

4A^ 



2e 



(sinh(Mfc) 



From Eq. (46) we see that at any k, the spectral function peaks at e* 



(46) 



(47) 



-e with the 



true quasiparticle peak corrected from the Fermi liquid value E^^i= ^k) as: = 
E[^ — 6^'^'= Z^Tk, at e* = — e""'-'. Simplifying, we find the EDC energy dispersion or spectrum 



El= Ck + eo-\^ [So + (1 - Z,) + Zin 



(48) 



We provide examples of this dispersion later in Fig. ([s]). 

For a given the magnitude of the spectral function at this peak is given by 



G*{k) 



n 



4Ao 



-Z. 



Phik) 



Zk 



(49) 
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The magnitude of 3?e G*{k) is a little smaller than the value arising at the FL solution 
e = 0. In Fig. ([T]) this is reflected in the line H2 that lies a little below HI [H] . 

The peak value p*Q{k) falls off with ^ 0, and is always smaller relative to the peak of 
the aux-FL peak value p*g{k). The ratio of the two peak values is given by 

P',(.k) 8A„ ^ ' 

We see below numerical examples of these functions. Fig. ^ illustrates the peak ratio and 
other features for a typical set of parameters. 



E. MDC or constant energy scans and energy dispersion. 

It is also useful to study the momentum distribution curves "MDC" obtained by fixing 
the energy x and scanning the energy [13]. In the model of a independent self energy, 



this is a particularly convenient strategy, and hence maximizing Eq. (37) at a fixed x we 
find the MDC energy dispersion or spectrum: 



f (x) =x-eo + Vr2(x) + (eo - 3?e<l>(x))2. (51) 



Thus ^*{x) is the peak position of in constant energy scans, whereas El in Eq. (48) 



represent peak position of energy at a fixed It is amusing to compare this with Eq. (48) 



Unlike Eq. (48), this formula is valid at all energies, not just near the chemical potential 
where the two agree closely. We will see below in Fig. ([s]) that this function is multivalued 
in a range of values of energy x leading to characteristic features of the spectrum. 



F. Numerical example of High Dimensional ECFL model: 

In this section we use a a rectangular band with height 1/{2W) and width 2W, and take 
W = .86 eV (i.e. lO'^K) as a typical value. In Ref. ^ a more realistic band structure is used 
as described in detail there. The model for the Fermi liquid introduced in Ref. ([1]) (Eq(24)) 
is given by the expression: 

T{x) = 7rp$(x) = nC^{x^ + r^je'^*^"'^"'^/"^ + r/, (52) 

with T = TiksT. We have added a scattering width rj as in Ref. ([3]), in order to account for 
scattering by off planar impurities. The real part of the self energy is found from the Hilbert 
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transform of p<i>{x), and is given by: 

^e<^(x) = C^n{x^ + r2)e-^*^"'+^'>/"^ Erfi(x/v^) - x e-^*"'/"% (53) 



where Erfi(x) = e*^ dt is the imaginary error function. A numerically small correction 

arising from rj is dropped for brevity. 

G. Typical parameters 

The same model is also used in the fit to experiments in Ref. (j3]) with a slight change of 
notation given by writing C$ — and ujc — )■ in terms of the high and low frequency 
cutoff frequencies loq and Qq. We use two sets of standard parameters 

Set-I: C$ = 1 eV-\ co^ = 0.25eV or loq = 0.5eV, fio = .318el^ 
Set-II: = 2.274 eV~\ uj^ = 0.568e\/ or coq = O.beV, fio = -1461/. (54) 

Set-I was used in Ref. ([T| for schematic plots employing a simple band density of states 
gsi^) = 2W®(^^ ~ ^^)' Set-II was used in Ref. employing a more elaborate dispersion 
described therein, to successfully fit data on various high temperature superconductors at 



optimal doping. The value of r] is displayed in different plots. In Eq. (39) the spectral 
function pg of the aux- PL is defined. The chemical potential is fixed by the number sum 
rule with ^ = e — fj, 

\ = j dx f{x) J de geie) Psi^- p,x), (55) 

where f{x) = (1 + e^^)~^ is the Fermi function. We now write the contributions from 
extreme correlations that are described in Ref. ([1]). The inelastic energy scale Aq is found 
from the sum rule: 

/oo r 
dx f{x) / de gsie) Pg(e - n,x) {e - p - x}. (56) 
oo J 

Thus at a given density and temperature n, T, the model has only two parameters Uc and 



C$ so that Aq is fixed from Eq. (56). We study the details of the spectra next. 
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V. THE SPECTRAL CHARACTERISTICS OF THE HIGH DIMENSIONAL 
ECFL MODEL EQ. ([37]) 



A. Global view of the spectral function 

We display in Fig. ([s]) the spectral function for the ECFL model Eq. (37) in 3-d plots and 
contour plots. Two distinct perspectives of the spectrum are found in the figure from the 3- 
dimensional plot and the contour plot. In both of these we see that the excitation are sharply 
defined only for a certain range near the Fermi energy, and then merge into the continuum. 
At higher binding energies, the spectrum again looks quite sharp. For comparison, in Fig. ([3| 
we also display the aux-FL spectral function. We note that the aux FL spectra also become 



sharp at higher binding energies. This sharpening is modeled by the Gaussian in Eq. (39), 
its basic origin is the decrease in the weight of physical processes capable of quasiparticle 
damping as we move towards the band bottom. In order to look more closely at the low 
energy part of the spectrum of the aux-Fl and the ECFL, we show in Fig. ^ the contour 
plots of both next over a smaller energy range. We see that viewed in this rather broad 
sense, dispersions of the aux FL and the model ECFl spectra share many characteristics, 
with somewhat different energy scales. However there are crucial differences that emerge 
when we look at the distribution of spectral weight that arises in the ECFL, where the 



caparison factor in Eq. (37) pushes weight to higher binding energies. This is reflected 
most significantly in the line shapes that we study below. Since we use the momentum 
independent self energy for the aux-FL in this model calculation, we obtain very detailed 
EDC and MDC plots below. However, it must be borne in mind that refined calculations 
within the ECFL framework must necessarily introduce some momentum dependence, and 
hence several details are likely to change, in particular the structure far from the chemical 
potential would change somewhat more. Our view is that this caveat apart, it is very useful 



to take the Eq. (37) seriously since it gives a simple framework to correlate different data. 



B. Dispersion relations in EDC and MDC 



In Fig. (|5| the EDC dispersion relation (i.e. locus of peaks of the spectral function at 
fixed ^, found by numerical maximization), is plotted versus ^ along with the MDC spectrum 



Eq. (51). We recall that the latter expression is exact at all ^,x, whereas Eq. (48) is not 
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quite exact for the EDC dispersion. For comparison we also show the corresponding figures 



for the aux- FL spectral function in Eq. (39), with the same parameters. The dispersion 
relations Eq. (48) is displayed in the inset of Fig. ([5]), where it is compared with the result 
of numerically maximizing the spectral function at a fixed C,- We see that Eq. (48) is only 
good for a range of energies near the Fermi energy. We see that both sets of spectra for the 
aux-FL as well as the ECFL model exhibit similar global features, but with different scales 
of energy. In both cases, the constant energy scans show a jump discontinuity, whereas the 
MDC spectra show a "S" type or re-entrant type behaviour. The origin of the latter is easy 
to see in the aux-FL, here a peak in — 3fJe$(a;) occurs at an energy approximately 2uq, so 
that as X decreases from zero, Caux-FL = ^ — 9fJe$(x) goes back up for a certain range. In 



case of the ECFL, Eq. (51) shows that for the energy scale enters the expression when 
r(a;) becomes comparable to — 3fte$(a;), and the net result is that the reentrant behaviour 
is pushed to lower binding energies. 



C. The energy shift 

The dispersion Eq. (48) corresponds to the ridge near the Fermi energy in Fig. ([3]). At 
low temperature, since F^^ ~ O(T^), the corrected quasiparticle energy is always less than 
E^^, so that there is always a leftward (i.e. red) shift of the dispersion, or from the hole 
(binding) energy point of view, we may say there is a blue shift. The peak shift is given by 



AE, = El - ^ = (1 - Z,)ik + £0 - V N + (1 - Zu) + ZlTl (57) 

which is a function of both T and k. Close to the Fermi energy this can be written as 

AE, = (58) 

At the Fermi momentum, this small shift is seen in Fig. ([T]) where the vertical line through 
X* is displaced to the left from the y axis. As long as A; ~ kp this shift is very small 



AE'fc ~ 0(T ), but as k moves away from kp the shift Eq. (57) grows with ^k- This departure 



makes the dispersion in Eq. (48) depart significantly from the bare dispersion ^ as we move 
away from kp- We see from Fig. ([s]) that the departure of the EDC peaks from the Fermi 
liquid are somewhat less pronounced than those of the MDC's, the latter is operationally 
called as the low energy kink. Our calculations therefore predicts the magnitude of the shift 
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Eq. (57) in terms of the energy scale Eq and the Fermi hquid parameter r^. This energy 
shift is therefore also a useful method for extracting the fundamental parameter Aq on using 
Eq. (|36|. 



D. Constant energy cuts or MDC line shapes 

We display the MDC line shapes in Fig. ([6]). Panel (A) shows the effect of the caparison 
factor ^1 — f + 'i^i^ " 3;)^, whereby the curves are skewed to the right, in contrast to the 
EDC curves that are skewed to the left. The latter important feature is also seen below in 
Fig. ([7]A) and noted in Ref. ([1]). Panel (B) shows the shallow peaks in the "S like" region 
of the energy dispersion seen in Fig. ([s]), and panel (C) shows the deep interior region where 
the peaks are more symmetric. 



E. Constant wave vector cuts or EDC line shapes 



The spectral function and the real part of the Greens function are calculated from Eq. (38 ) 



and Eq. (37). We display the EDC line shapes in Fig. ([T]). Panel (A) gives an overview of the 
spectral shapes for wave vectors near the Fermi surface, displaying a left skewed peak that 
falls rapidly in intensity as it broadens. This behaviour is of great interest since it captures 
the experimental features in high Tc systems, as elaborated in Ref. Panel (B) shows 
the spectra at higher binding energies, where a feature at lower energies begins to disperse 
significantly with C,- It is evident that these two sets of dispersing features correspond to 
the two branches that are seen in the 3-d plots and contour plots Fig. ([3]). The inset in (B) 
shows the behaviour of the aux- FL, where the two features are again seen but with different 
rates of intensity change. 

We now turn to the task of understanding the reconstructed Dyson self energy that leads 
to the above electron spectral functions. In Fig. ([s]) we show the spectral function pc at 
various values of the energy ^ at T = 300K. The Fermi liquid spectra at the same values 
of parameters are also shown for comparison. The Dyson self energy Pt,{C,,x) necessary to 



produce these spectral functions is found using Eq. (23) and displayed in Fig. ([8j) at two 
temperatures. The object pj:{C,,x) has a distinctive minimum for each ^ that shifts to the 
left along with the energy ^ that tracks the peaks in the physical spectral function Pg{^, x) 
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from Eq. (23). It also shows the asymmetry between energies above and below the chemical 
potential that we noted at ^ = in Fig. At the Fermi energy Py:{^,x) is displayed in 
Fig. ([9]) over a large scale. 



F. The reduced line shape function 



An interesting aspect of the ECFL model Greens function Eq. (37) is the change in 



shape of the peaks as we leave the Fermi surface, so that the quasiparticles become hard to 
define at some point. This change in shape can be formulated neatly in terms of a single 



dimensionless parameter Qk that we now define and explore. We examine Eq. (46) around 
its peak by writing 

e = e* + cosh(ufc)e, (59) 
so that pQ^"'^{ik^(^) = PG'(^)7(*5fc5 e), with a characteristic line shape function 7 given by 

Qk{l - e) 



7(Qfe,e) 



Qfc(l-e) + e2 



with 



Qk 



(60) 



(61) 



cosh (uk) 

The parameter Qk goes to zero near the Fermi surface at low T since Wfc — )■ 00, but at higher 
binding energies increases Qfc — 2. 

As we get deeper into the occupied states ^ <^ 0, we find a remarkable change in shape 



of the spectral functions. This is illustrated in Fig. (10) where we plot 7 of Eq. (60) after 



normalizing to unit area. In order to have a well defined quasiparticle type peak in pc for 
e ~ 0(1), Qk must be small enough. This translates to the requirement of ^k being close to 
the Fermi surface. By setting Q ~ 1 as the condition for losing a peak in the spectrum we 
obtain the condition 



£0 + (1 - Zk)^k 



ZkT 



fcJ- k, 



(62) 



beyond which it is meaningless to talk of quasiparticles. This gives £0 as a rough character- 
istic scale for the disappearance of the quasiparticle peaks. 



Fig. (11) illustrates the change in shape somewhat differently by normalizing all curves 



to unity at the peak as in Eq. (60). The peak at e = is sharp and quite symmetric for 



Q <^ 1, and becomes broader and more left skewed as Qk increases towards its maximum 
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value 2. Attaining the maximum value is possible in principle, requirement being = or 



from Eq. (44): 

,2 



^ 5Re (63) 



4 - 2n 

Unless Aq is very small, this condition is hard to satisfy. If this possibility is achieved then 



there are several interesting consequences. Firstly we note that from Eq. (49), the magnitude 
of the spectral function at = becomes insensitive to disorder and temperature etc. Its 
magnitude n^Zfc/(87rAo) that should be useful for finding Aq. If this is approximately 
satisfied, then the peak structure loses meaning and the spectrum is essentially fiat. Taking 
e = — 1, the fall off from the peak value of unity is 80%, and the spectrum becomes essentially 
featureless. 



G. Skewness parameter of the spectrum 



We now estimate the skewness of the spectrum. The function Eq. (60) drops to half its 



peak value at two values of eL,R to the left and right of the peak given by: 



eL,R = ~\QkT\jQk + lQl (64) 

For small Q ^ 1, the (dimensionless) width of the peak is small, it increases with Q as 
discussed further below. We define a dimensionless skew parameter in terms of the energies 



e*,eji^L (rather than the e's that are related via Eq. (59)) 



<ik) = , = tanh(Mfc) - 1 + v/(2-tanh(Mfc))2-l, (65) 

with the property that near the Fermi level when Uk — )■ +oo the variable k — )■ 0, and we 
get a symmetric curve about the maximum. On the other hand for deeper occupation Uk 
decreases in magnitude towards zero, driving k, — 0.732 and gives a curve that is increasingly 
biased to the left. The asymmetry k grows as 0{T'^) at low temperatures, and it is rather 
large at room temperature. As a rough estimate the quasiparticle peak is lost when Qk ~ .5 
where Uk ~ .98 and K,{^k) ~ -5. This loss of quasiparticle peak structure, skew factor and 
its experimental signature is studied in greater detail in Ref. See Fig. ^ for a typical 
plot of skew parameter k and Q versus the hole binding energy E^. 
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VI. SELF ENERGY OF THE SELF ENERGY AND A MORI TYPE PROCE- 
DURE. 



Since the construction given above generates ^{z) from G{z) given only the representation 
Eq. 0, we can as well repeat the trick. Since J]{z) satisfies Eq. ( |18| ) with a density pj:{x) 
that is assumed known, and is analytic in the complex z plane with al/z fall off, it satisfies 



the necessary conditions for a further decomposition. Consulting Eq. (23), we write down 
by inspection 



z-Ej:-EW{z) 

eW(.) = [dx (66) 
J z-x 

The constants as = j dx ps{x) and Ej] = j dx xpj:,{x)/aj] are known through p-£{x). They 
may in practice be conveniently determined in terms of the moments of the Greens function 
[3-[9] in applications. The spectral function is given by 



Qs PT.{X) 

(7r2ps(a;))2 + (3?e E(x))' 



Comparing this representation with Eq. ( plj ) we note the formal similarity between Pg{x) 
and p-£(i){x). Thus for a Fermi liquid with momentum independent self energy, its next self 
energy resembles closely the spectral function pc especially at the Fermi energy. 

We follow up briefiy on the above amusing observation, and obtain a hierarchy of self 



energies starting from an initial self energy given by the spectral representation Eq. (18) 



and Eq. (66). This process parallels the continued fraction representation of analytic func- 
tions and seems intimately related to the formalism developed by Mori[71 18]. The latter is 
expressed in the language of projection operators for Liouville operators that is less straight- 
forward than our simple treatment. 

In order to conform to the notation popular in the Mori formalism, we will express the 
variables in Laplace representation rather than the one used above with complex frequencies. 
Let us consider the thermal and temporal correlation function for two operators A and B 
in Schrodinger time and its Laplace transform: 

/•/3 POO 

CAB{t) = / dr (A(t-zr)5(0)), Cab{s)= e''' CAsit) dt (68) 
Jo Jo 
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In the standard case we find A = B\ where the product is real and also positive [16]. We 
see that the Laplace transform function satisfies an integral representation 

Cab{s) = [ dv PAeLJ.^ with a real density given by 
./-^ s -iv 



Pab{^) = y^ ^" ^" (ra|y4|m) {m\B\n) 6{en - £m - i^)- (69) 



This object is closely connected with the correlation functions used in Eq. ([T]) and Eq. ^ 
by using the fluctuation dissipation theorem. Following Mori we write down a relaxation 
function with the normalization property Zo{0) = 1, and its Laplace transform 

= c^y ^°<^' = y_ J^u (™) 

Here the real density Pq{v) = ^^pab{.^), satisfies the normalization condition 
du Pq^u) = 1. Using the identity 

7T 0{U — Vj — iV- 



with V denoting the principle value, an inverse relation expressing Pq{i') = Zq{0^ + 
follows. In order to find a Dyson type representation for po, following Eq. ([s]) and Eq. ([9]) 



we take the "equation of motion" by multiplying Eq. (70) by s and write 



oo 



sZ,{s) = l + i I du^'-^^Y.is) Z,{s). (71) 



oo 



The Dyson form of self energy now emerges and we obtain: 



1 r dv 

Zo(s) = with Y^is) = (72) 

As s — 7- oo, the function Yq{s) tends to Wi, with a real frequency Ui given by 

Ui = dv v po{v). (73) 



Hence the function (Yo{s) — Ui) falls off as ^ as s — )■ oo. It is analytic everywhere except on 
the imaginary s axis. It therefore has a representation 

YJs) - ui = tai r dv (74) 
y_oo s-tv 

with a real density aipi{v) = ^Qm{YQ{0'^ + iv)}. With this we may write 



s — i ui + ai J dv 
23 



The real number ai is found using the convention that pi(i^) is normahzed to unity. We 



may express pi solely in terms of the lower density poii^) by using Eq. (70) as 



Po[u) 



«iPi(w) = ^. (76) 

n'pl{u) + {H[po]{u)}' 



We determine ai from Eq. (76) by integrating over u and using the unit normalization of 
Pi{u). It is evident from Eq. (76) that for the physically important case of a real and positive 
initial density Po('^), the derived density piiy) is also real positive. 

This scheme is clearly generalizable to higher orders, and we simply iterate the above 
process. The answers may be written down by inspection as follows. 

Zj{s) = dv — — — , with normalization: / dv Pj{v) = 1. (77) 

J -oo S — IV J -oo 

These satisfy a recursion relation: 

ZAs) = : ^— , (78) 

where 

/oo 
du u pj{u), (79) 
-oo 

and ctj+i as well as Pj+i{i^) are defined through 

p,+i (u) = /^'^^^ . (80) 

Note that the numbers aj as well as uj are real, and for all j, the densities Pj(z^) are positive 
provided the the initial density Po('^) is positive. This situation arises when the initial 
operators B = A\ as mentioned above. 



It is clear that Eq. (76) is the precise analog of the relation Eq. (23) for the Greens 



function. The hierarchy of equations consisting of Eq. (45-49) constitutes an iteration scheme 



that starts with j = correlation function in Eq. (70). This is a forward hierarchy in the 
sense that successive densities at level j + 1 are expressed explicitly in terms of the earlier 
ones at level j. In the reverse direction it is rather simpler since level j is explicitly given 



in terms of level j ' + 1 by Eq. (78). The use of this set of equations requires some a priori 
knowledge of the behaviour of higher order self energies to deduce the lower ones. Standard 
approximations [7J consist of either truncation of the series or making a physical assumption 
such as a Gaussian behaviour at some level and then working out the lower level objects. 
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Our object in presenting the above working is merely to point out that this iterative scheme 
is in essence a rather simple application of the self energy concept described above, with the 
repeated use of Eq. (23). 



VII. SUMMARY AND CONCLUSIONS 

A new form of the electronic Greens function, departing widely from the Dyson form 
arises in the extreme correlation theory of the t-J model. Motivated by its considerable 
success in explaining ARPES data of optimally doped cuprate superconductors |3], we have 
presented in this paper results on the detailed structure of this Greens function and its 
spectral function. An illustrative example is provided, complete with numerical results, 
so that the novel line shape and its dependence on parameters is revealed. We have also 
presented a set of explicit results on the Mori form of the self energy that holds promise in 
several contexts. 

VIII. ACKNOWLEDGEMENTS 

This work was supported by DOE under Grant No. FG02-06ER46319. I thank G-H. 
Gweon for helpful comments and valuable discussions. I thank A. Dhar and P. Wolfle for 
useful comments regarding the Mori formalism. 



[1] B. S. Shastry, [ arXiv: 1102.28581 (2011), Extremely Correlated Fermi Liquids Phys. Rev. Lett. 

107, 056403 (2011). 
[2] B. S. Shastry, Phys. Rev. B 81, 045121 (2010). 

[3] G.-H. Gweon, B. S. Shastry and G. D. Gu, Extremely Correlated Fermi Liquid Description of 
Normal State ARPES in Cuprates, arXiv:1104. (2011), Phys. Rev. Lett. 107, 056404 (2011). 
[4] D. E. Logan, M. P. Eastwood and M, A. Tusch, J. Phys. Cond. Mat. 10, 2673 (1998). 
[5] Z. Wang, Y. Bang, and G. Kotliar, Phys. Rev. Lett. 67, 2733 (1991). 

[6] J. M. Luttmger and J. C. Ward, Phys. Rev 118, 1417 (1960), J. M . Luttinger, Phys. Rev. 

119, 1153 (1960); Phys. Rev. 121, 942 (1961). 
[7] H. Mori, Prog. Theor. Phys. 33, 423 (1965); Prog. Theor. Phys. 34, 399 (1965). 

25 



[8] M. Dupuis, Prog. Theor. Phys. 37, 502 (1967) 

[9] J. J. Deisz, D. W. Hess and J. W. Serene, Phys. Rev. 55 2089 (1997). 

[10] A. A. Abrikosov, L. Gorkov and I. Dzyaloshinski, Methods of Quantum Field Theory in Sta- 
tistical Physics , Prentice-Hall, Englewood Cliffs, NJ (1963). Our poiCk^x) corresponds to the 
combination of spectral functions A{k, x) and B(k, x) used here. 

[11] The positive part constraint on the right of the following equation can be often omitted, we 
found that it is violated very slightly ~ 3% in many cases. We omit it for simplicity in the 
following. 

[12] Shifting E hy a. constant is also possible but the optimal choice is the one made here. 

[13] The notation is simplified from that in Ref. (1) by calling the extremely correlated Greens 
function as G rather than Q and the overbar in the self energy $ is omitted. 

[14] At T = this finite value persists and 3?e G^'^°'^{^k: x*{^k)) does not diverge, so that one might 
be concerned that the Luttinger Ward volume theorem is being disobeyed. In comparison note 
that the standard FL Greens function behaves in a slightly different way, at any finite T and 
^k, at the energy x = E^^ we find both a peak in the spectral function pg{^k, E^^) and a 
zero of the Ke g{^}^,E^^), whereas at T = 0, we find a delta peak in the spectral function 
Pgi^k, E^^) and a pole of the Ke g{^k, E^^). However we see that the FL divergence of the 
real part of G of a typical Fermi liquid does occur, but displaced by a very small energy scale 
O(T^). The peak positions are displayed in Fig. ([l]) at a high enough temperature so that the 
features are distinguishable. 

[15] The process of scanning used here differs slightly from the true MDC's, where one scans 
the wave vector k rather than the energy ^fc, but is more convenient here. 



[16] The spectral function PAA^i-'^) defined below in Eq. (69), is positive as well. However that 
condition can be relaxed and we can do with less provided that the spectral function in 



Eq. (69) is real (rather than positive). For this to happen, we may allow for A B"^ , but 
in this case assume that both the matrix elements (n|^|m) and (m|i3|n) are real numbers or 
imaginary numbers so that the product is real. These conditions correspond to both operators 
operators A,B being Hermitean (or anti Hermitean), and the absence of magnetic fields so 
that the wave functions may be chosen to be real. Thus we will assume the reality of the 
product of the matrix elements. 



26 




-3 






FIG. 3: The ECFL (top left) and the auxihary FL spectral function (bottom left) at density 
n = .85, T = 180K, Aq = 0.12 eV, ij = 0.12 eV, and the other parameters are from Set(II) in 



Eq. (54). Here ^ and x are in units of eV. In the ECFL curve on left, it is seen that the excitations 
near the Fermi energy become broad and dissolve into the continuum at an energy ~ —.2 eV, 
and reappear as sharp modes at a deeper binding energy. In the auxiliary FL, the excitations 
near the Fermi energy remain sharp and extend to lower energies than in the ECFL curves. The 
contour plots of the same functions in the right panel (top ECFL and bottom auxiliary FL) gives 
a complementary perspective of the spectrum. The two superimposed solid lines at top right are 
from Curves (I) and (II) of Fig. ([s]) and at bottom right Curves (III) and (IV) of Fig. ([s]). 
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-0.7-0.6-0.5-0.4-0.3-0.2-0.1 0.0 -0.7-0.6-0.5-0.4-0.3-0.2-0.1 0.0 

FIG. 4: The contour plots of the aux FL (left) and the ECFL model from Fig. ([3| (right) with the 
same parameters as in Fig. @ but over a smaller energy window. We superimpose the constant 
wave vector dispersion and MDC dispersion, with a value of = .12 common to the contour plots. 
The energy scale of the feature near the chemical potential is considerably reduced in the ECFL, 
and the "jump" in the EDC dispersion occurs at roughly half the corresponding energy in the 
aux-FL. 
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FIG. 5: Energy dispersion curves in the ECFL and the aux FL models. Here the parameters are 



from Set(II) in Eq. (54), with n = 0.85, T = 180K. With r/ = 0.12, curves I and II are the peaks 



in constant wave vector and constant energy scans of the spectral function Eq. (37), and curves 



III, IV are corresponding figures for the aux FL in Eq. (39). The inset compares Eq. (48) (the 



truncated curve) with the exact locus found by numerical maximization. 
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FIG. 6: MDC line shapes at different values of energy x displayed in each curve. Here the 



parameters are from Set(II) in Eq. (54), with n = 0.85, T = 180K and rj = 0.12. Panel (A) 
corresponds to x close to the chemical potential. It is interesting to note that the curves are 
skewed to the right, thus mirror imaging the leftward skew seen in the constant ^ (EDO) scans 
below Fig. ([7]A), in a comparable range of energy and wave vector. Panel (B) corresponds to the 
mid energy range, within the re-entrant range of x from Fig. ([s]) or Fig. ([s]), with the counterintuitive 
movement of the shallow peak to the right with increasing x. Panel (C) corresponds to the second 
set of maxima in Fig. ([s]) far from the chemical potential, where the curves are quite symmetric. 
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FIG. 7: EDC line shapes at different values of energy ^ displayed in each curve. Here the 



parameters are from Set(II) in Eq. (54), with n = 0.85, T = 180K and i] = 0.12. Panel (A) 
corresponds to ^ close to the chemical potential. Note that the curves are skewed to the left, 
i.e. a mirror image of the rightward skew seen in the constant x MDC scans above Fig. ([g]), in a 
comparable range of energy and wave vector. Panel (B) corresponds to the higher energy range, 
and we see that only one broad maximum is found at a given ^. The inset of (B) shows the aux-FL 
constant ^ scans for the same range, here each ^ results in a pair of maxima, originating from the 



functional form of the self energy in Eq. (37). 
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FIG. 8: Top Left panel The density n = .85, temperature T = 300K, Aq = .0678 and other 



parameters are from Set(I) in Eq. (54). From left to right pg{^,x) for energies in units of eV: 
^ = -0.1,-0.075, -0.05, -0.025,0., 0.025, 0.05. Top Right panel The spectral function pg{^,x) 



from Eq. (39) corresponding to the same ^ as in the left panel. The difference in the line shapes 



becomes clear when we examine the Dyson self energy that produces these curves. Bottom left 
panel The panel shows the spectral function at T = 300K for the inferred Dyson self energy pt,{x) 



from Eq. (23) and Eq. (38), Eq. (37) for the same energies. The dashed line is the input Fermi 



liquid spectral function p^{uj) at the same temperature. Bottom Right Panel Temperature 
T = 150, Aq = .0642 and the identical data as in the bottom left panel. 
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FIG. 9: With n = 0.85 and T = 300i^, and other parameters are from Set(I) in Eq. (54). The 
spectral function for the inferred Dyson self energy /Os(C = 0,x) using Eq. ( p3| ) and the the input 
Fermi liquid spectral function p^{x) over a larger energy range. Note the distinctive asymmetry 
in shape of below above the Fermi energy. 
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FIG. 10: The spectral shapes possible are seen by plotting the shape function different values of 



the parameter Q. In this curve 7 is the 7(Q, e) of Eq. (60) normalized to unit area in the natural 
interval [—1,1] for the variable e. 
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FIG. 11: The same shapes as in Fig. (10) but now 7 is normahzed to unity at the peak as in 



Eq. (60). The sharp peaks for small Q < .25 flatten out as Q increases with a left skew asymmetry 



that is characteristic of this functional dependence. 
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